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Isothermal crystallization of the mono-hydroxyl alcohol n-butanol was studied with dielectric 
spectroscopy in real time. The crystallization was carried out using two different sample cells at 15 
temperatures between 120 K and 134 K. For all temperatures, a shift in relaxation times to shorter 
times was observed during the crystallization process, which is characterized by a drop in relaxation 
strength. The two different sample environments induced quite different crystallization behaviors, 
consistent and reproducible over all studied temperatures. An explanation for the difference was 
proposed on the background of an Avrami and a Maxwell-Wagner analysis. Both types analysis 
suggest that the morphology of the crystal growth changes at a point during the crystallization. 
The differences between the cells can be explained by this transition taking place at different times 
for the two cells. 


All liquids can be supercooled [1-3]. In fact, crystal¬ 
lization rarely takes place exactly at the melting temper¬ 
ature upon cooling, because the crystal nuclei formed dis¬ 
solve before they grow to a stable size [4]. Some liquids, 
like water, crystallize readily at moderate supercooling 
and need fast quenching below the melting temperature 
in order to avoid crystallization. Others, like the proto¬ 
type glass-former glycerol, supercool easily and require a 
careful protocol to crystallize [5]. But the true thermo¬ 
dynamic equilibrium state for all supercooled liquids and 
glasses is unarguably the crystal, and thus crystallization 
is their eventual inevitable fate. 

For many applications the life-time of the glassy or 
meta-stable liquid state is a key issue and the ability 
to predict and control crystallization properties is de¬ 
sired. But understanding the crystallization process is 
also interesting from a fundamental point of view. Crys¬ 
tallization studies are however difficult to carry out in a 
controlled and reproducible manner, because many fac¬ 
tors influence the initiation and course of crystallization, 
such as sample preparation, thermal history, presence of 
impurities, container geometry, etc [5-7]. 

Dielectric spectroscopy is a convenient and fairly com¬ 
mon probe for studying crystallization kinetics in real 
time, see e.g. Refs. [8-13]. The sign of crystallization 
is a decrease in intensity of the signal, and the crystal 
fraction/concentration of the sample is often obtained 
by assuming that the relaxation strength is proportional 
to the volume fraction of liquid in the sample [14-16]. 
But the microscopic interpretation of dielectric spectra 
is not straight-forward [17], especially when studying a 
heterogeneous mixture. 

We present a thorough dielectric study of the tem¬ 
perature dependence of the crystallization process in 
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n-butanol. Supercooled n-butanol has an intense low- 
frequency dielectric signal - the so-called Debye-process 
- characteristic of many monohydroxyl alcohols [18]. It 
also exhibits a slow crystallization process upon reheat¬ 
ing after a rapid quench below Tg [19, 20]. At 10 K 
above the glass transition temperature the crystalliza¬ 
tion can take several days to finish, and at temperatures 
close to the glass-transition, the crystallization process 
stops before the sample is fully crystallized [21]. These 
two properties, a large dielectric signal and slow crys¬ 
tallization, makes n-butanol an ideal candidate for moni¬ 
toring isothermal crystallization in real time by dielectric 
spectroscopy. Crystallization of n-butanol has previously 
been studied with x-ray diffraction [19, 20], infrared spec¬ 
troscopy [21, 22], calorimetric methods [23], and phase 
contrast microscopy [24]. However, no systematic study 
of the temperature dependence and reproducibility of the 
crystallization kinetics has been reported. We used two 
different dielectric measuring cells with different geom¬ 
etry and different electrode material in an attempt to 
disentangle effects that are intrinsic to the sample and 
effects that are due to macroscopic boundary conditions. 

I. EXPERIMENT AND MATERIALS 

All measurements were carried out in the same exper¬ 
imental set-up (described in detail in Ref. 25), includ¬ 
ing a custom-built nitrogen cryostat capable of keeping 
the temperature stable within 50 mK over weeks. Two 
different sample cells were used. Cell A is a 22-layered 
gold-plated parallel plate capacitor with 0.2 mm between 
each set of plates and a geometrical capacitance of 64 pF. 
Each plate is a semi-circle which can be rotated to over¬ 
lap each other (identical to the capacitors used in old 
radios). Cell B is a parallel plate capacitor with circu¬ 
lar beryllium-copper plates separated by 50 jim. sapphire 
spacers. The cells are sketched in Fig. 1. 
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FIG. 1: Schematic drawings of the two cells used for the 
crystallization studies. Cell A is a multilayer capacitor 
(variable capacitor type used in old fashioned radios) 
with a geometric capacitance of 65 pF. Cell B is a 
two-plate capacitor with beryllium copper electrodes 
and 50 /im sapphire spacers. Cell B has a geometric 
capacitance of 17.6 pF. 


The sample n-butanol (T^ = 183i^, Tg « llOA") was 
purchased from Sigma Aldrich’s at > 99.9% purity and 
used without further purification. 

The same protocol was followed for each crystalliza¬ 
tion measurement; a new sample was quenched to 85 K 
(roughly 25 K below Tg) and kept at this temperature for 
(at least) 120 minutes, then heated to the target temper¬ 
ature where the crystallization process was followed. The 
heating took less than five minutes. Frequency scans were 
made continuously as soon as heating from 85 K initiated 
and until no further changes in the spectrum occurred. 
The frequency range of the scans was adjusted for each 
temperature to keep the scan as short - and thereby as 
fast - as possible, while still keeping both Coo and in 
the frequency window. The sample cells were emptied 
and cleaned between each measurement. 

Isothermal crystallization was followed at 15 different 
temperatures between 134 and 120 K. The exact tem¬ 
peratures in Kelvin are: 134, 133, 132, 131, 130, 129, 
128.5, 127.5, 126, 125, 124, 123, 122.5, 121, 120. The 
133 K measurement has only been done with cell A, while 
the 132 K measurement has only been done with cell B. 
For reference, a fully crystallized sample was made by 
quenching to 85 K and reheating to 170 K, and a spec¬ 
trum was measured at all temperatures included in the 
study. 


II. EXPERIMENTAL RESULTS AND DATA 
ANALYSIS 

A. Phenomenological fits of the spectra 

In n-butanol there are three visible processes - Debye, 
alpha, and beta process - in the measured frequency win¬ 
dow, see Fig. 2(a). The crystallization, signaled by a de¬ 
crease of relaxation strength, also induces a shift in the 
loss peak for the three processes. To quantify how the 
crystallization influences each of these processes, we fit¬ 
ted the spectra to a sum of three relaxation processes. 
Since the processes are not well separated, we aimed at 
limiting the number of free fitting parameters by the fol¬ 
lowing procedure: First the Debye process is fitted by 
a Cole-Cole function. The Debye process broadens dur¬ 
ing the crystallization and thus a pure exponential func¬ 
tion would not give a good fit. The result of the Debye 
fit is then subtracted from the data and the alpha and 
beta processes are fitted simultaneously as a sum. The 
beta process is fitted to a Cole-Cole function with a fixed 
shape parameter, (3 = 0.45. The alpha process is fitted to 
a dielectric version of the Extended Bell (EB) model (see 
Ref. 26), in which the imaginary part of the dielectric 
constant is given by: 

^EB = I (1) 

1 

{l + iu)Tay +fca(zwra)“" 

where Ae is the relaxation strength and is the relax¬ 
ation time, ka controls the width of the peak, and a gives 
the high-frequency power law behavior of the alpha peak. 
This model for the alpha relaxation gives good fits even 
when the shape parameters are fixed such that only the 
relaxation strength and relaxation time is fitted. The 
slope parameter was fixed to a = 0.5 [27] and = I- 
The Cole-Davidson function resulted in poorer fits, even 
with the shape parameter varying freely. Thus, we fit¬ 
ted the imaginary part to a total of seven parameters: a 
relaxation time, r^,, and a relaxation strength, Ae^,, for 
each process {x denoting Debye, alpha, or beta), and a 
broadening parameter for the Debye process. 

The procedure gives excellent fits of over the frequency 
range explored as shown in Eig. 2(a), where the first 
dielectric spectrum at every annealing temperature is 
shown as well as the fits resulting from the described 
fitting procedure. The inset shows a spectrum at 130 K 
with each of the individual fitted relaxation processes. 

The same fitting procedure was applied to the isother¬ 
mal crystallization spectra to study the temporal evolu¬ 
tion of the three processes during crystallization. This is 
shown in Eig. 2(b) at 130 K using cell B. As the crys¬ 
tallization proceeds the strength of the Debye and alpha 
relaxation processes decrease continuously to disappear 
entirely by the termination of the crystallization process. 
This is also to be expected, since there should be no 
large-scale rearrangement of the molecules in the crystal. 
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FIG. 2: Fitted spectra and time and temperature 
dependence of the fitting parameters, (a) The first 
measurement at every temperature, presumably before 
crystallization initiates. The magenta lines are fits 
described in the text. The inset show the three 
relaxation processes constituting the fit: the Debye 
(green), alpha (red) and beta (blue) processes, (b) 
Selected spectra during crystallization at 130 K, using 
cell B. Again, magenta lines are fits. The inset show the 
dielectric loss of the last measurement (blue), together 
with the measurement of the full crystal (black), on a 
linear scale. The ’spike’ at 100 Hz is an experimental 
artifact deriving from a non-perfect match between two 
voltmeters, (c) The fitted relaxation times together 
with the crystallization time. The lines are linear fits, 
(d) Parameterised plot of the fitted relaxation 
relaxation frequency, / = 1/r, against the fitted 
relaxation strength Ae, for each of the three processes. 
The two cells display quite different behavior during 
crystallization. 


The beta relaxation process, however, remains active by 
the termination of the crystallization process. The inset 
of Fig. 2(b) shows the last scan at 130 K together with 
a measurement of the full crystal at the same tempera¬ 
ture, clearly demonstrating that the crystallization pro¬ 
cess stops before the sample is fully crystallized and that 
there is still some molecular mobility left. The results 
from the fitting routine establish the general behavior 
of the spectra during the crystallization, but we refrain 
from analyzing the finer details, especially towards the 
end of the crystallization process, where Debye and al¬ 
pha process has vanished and the fits become unreliable 
(the lowest curve in Fig. 2(b)). 

Using the relaxation strength, Ae, as an indicator of 
the degree of crystallinity, we define a characteristic crys¬ 
tallization time as the time the for Ae to decay to half 
of its initial value. Figure 2(c) shows the crystalliza¬ 
tion time derived from both Debye and alpha relaxation 
strength as a function of inverse temperature, and it is 
evident that the two measures are identical within the 
accuracy of our measurements. 

Along with the crystallization time, we show the re¬ 
laxation times obtained from the fits to the first (un¬ 
crystallized) spectrum. Clearly, the Debye, alpha, and 
beta processes as well as the crystallization process are 
all slowed down with decreasing temperature and con¬ 
sequently the different characteristic time scales would 
all appear to be correlated (at least in this temperature 
range), but it does not necessarily imply any causation. 
For the studied temperature span, all the shown time 
scales are Arrhenius within the noise, although with very 
different pre-factors. 

Focusing now on how the two sample environments in¬ 
fluence the crystallization process, we define the relax¬ 
ation frequency as the inverse of the fitted relaxation 
times, / = 1/r. The evolution of both the relaxation 
strength and relaxation frequency differ for the two cells. 
Fig. 2(d) shows a parameterised plot of the fitted re¬ 
laxation frequency (normalized to the initial value) as a 
function of the fitted normalized relaxation strength for 
each of the three processes from measurements at 127.5 K 
with both cell A and B. In cell A, there is a shift in the 
relaxation frequency, /, quite early in the crystallization 
process, then it remains relatively unchanged for all three 
relaxation processes until an increase sets in again to¬ 
wards the end of the crystallization. In cell B, the shift 
in relaxation frequency is more gradual. For the Debye 
process the shift is monotonous, but for the alpha and 
the beta process the shift displays a non-monotonous be¬ 
havior. The curve peaks in Fig. 2(d) occurs roughly the 
same waiting time for the alpha and beta process. 

The different progress of the crystallization process for 
the two cells suggests a macroscopic/mesoscopic rather 
than microscopic explanation since a slight difference in 
sample geometries is not expected to affect the behavior 
of individual molecules. 

The full set of fitted parameters normalized to the ini¬ 
tial value is shown as a function of waiting time in Fig. 3. 

















4 





Normalized Ae 


FIG. 3: The fitted relaxation strengths and relaxation times for isothermal crystallization curves og n-butanol for 
temperatures between 120 K and 134 K (see text for details on the fitting procedure). The colours of the curves 
indicate the temperatures with blue being the lowest temperature and red the highest. 133 K has only been done 
with cell A while 132 K has only been done with cell B. Notice that the crystallization process takes longer to finish 
in cell B.(a+d) The normalized relaxation strengths defined as Ae(t)/Ae(0), of the Debye, alpha, and beta processes 
for cell A (a) and B (d) as a function of logarithmic time, (b+c+e+f) The logarithm of the normalized relaxation 
frequency (log [/(t)//(0)], with / = 1/r ) as a function of the normalized relaxation strength for cell A (b+c) and 

cell B (e+f). 


The colors of the curves indicate the temperature with 
blue being the lowest (120 K) and red being the highest 
(134 K). In both cells, lower temperatures lead to longer 
crystallization times, as was shown in Fig. 2(c). 

For cell A (Fig. 3(a)) the fitted relaxation strengths 
for each of the three processes appear similar, except at 
long waiting times, where the beta relaxation strength 
levels off at ^ 0.08 instead of decaying all the way to 
zero. The final level for the beta relaxation is marked 
by a dashed line in the third panel of Fig. 3(a). For cell 
B (Fig. 3(d)), the Debye and alpha relaxation strengths 
follow each other until roughly halfway through the crys¬ 
tallization, where a shoulder emerges in the alpha relax¬ 
ation strength curve, which then proceeds like a two-step 
relaxation. As in cell A, the beta relaxation strength 
does not decay to zero and levels off at the same value as 
for cell A. Comparing Fig. 3(a) and (d), we see that the 


curves for cell B are significantly more stretched than the 
corresponding curves for cell B, which means that crys¬ 
tallization proceeds at a consistently slower rate in cell B 
compared to cell A. Consequently, our definition of crys¬ 
tallization time may give roughly the same for the two 
cells, however total crystallization time is much longer in 
cell B. 

Figure 3(b+c) and (e+f) show the parameterized plot 
introduced in Fig. 2(d) for all temperatures, i.e., nor¬ 
malized relaxation frequencies as a function of relaxation 
strength, but now separating the two cells and the Debye 
and alpha process. (The relaxation frequency of the beta 
process is not shown, because it does not vary in a sys¬ 
tematic way, making further interpretation unjustified.) 

In cell A, the general behavior is that the relaxation 
frequencies have a slight shift to higher frequencies at 
the onset of the crystallization, but only increase a little 
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during the remainder of the crystallization. In cell B, the 
relaxation frequency of Debye and alpha processes do not 
change in the beginning of the crystallization process, but 
shifts gradually to higher frequencies. A ’bump’ occurs in 
alpha relaxation strength around a normalized relaxation 
strength of 0.5, showing that the behavior observed in 
Fig. 2(d) is general. 

The differences between the two cells are thus repro¬ 
duced for all the studied temperatures. One mesoscopic 
explanation for the observed difference in the evolution of 
the relaxation strength and relaxation times for the two 
cells could be that the two cells induce different kinds of 
crystal growth. 

B. Maxwell-Wagner analysis 

For heterogeneous material, a difference in the con¬ 
ductivity of the different domains in the material leads 
to build-up of charges at the interfaces between domains. 
This gives rise to a polarization effect known as Maxwell- 
Wagner (MW) polarization [17, 28]. 

In the present case the heterogeneity is caused by the 
formation of crystallites in the sample. As domains of 
crystal grow in the liquid the dielectric constant for the 
composite will change. The details of the change will 
depend on the difference between the dielectric constant 
of the liquid, e/, and that of the crystal, ec, the shape of 
the crystal domain, and the volume fraction taken up by 
the crystal. 

The two simplest cases of crystal domains growing in 
the liquid are that of a crystal layer growing from one 
(or both) of the electrodes, a heterogeneous nucleation 
picture, and the case of crystal spheres in a liquid matrix, 
a homogeneous nucleation picture. 

In the first case, no approximation is involved in de¬ 
riving the expression for the composite dielectric con¬ 
stant. The two materials (liquid and crystal) in a layered 
construction are simply modeled by two capacitors con¬ 
nected in series. Thus the resulting composite dielectric 
constant is given by 


(l-<^slab)ec + <^slabe; ’ '' 

where 0siab = d^jd is the relative thickness of the crystal 
layer. Since ei and €c can be measured independently 
this model has one free parameter (assuming the distance 
between the electrodes is fixed, or equivalently that total 
thickness of crystal and liquid layer is unchanged during 
crystallization). 

Inserting the measured spectrum of ei (at time t = 0 
before crystallization initiates) and Cc (spectrum of the 
fully crystallized sample), this model produces a large 
shift in the peak frequency for even small values of 0siab- 
Thus, the model is unable to capture both the observed 
decrease in relaxation strength and the shift in peak fre¬ 
quency at the same time, and this scenario alone is not 
sufficient to explain what we observe. 


In the case of crystal domains dispersed in a liquid a 
mean-field approximation is used to arrive at the com¬ 
posite dielectric constant [29] 

_ 2e/ + ec — 20(e/ — Cc) . . 

^comp — ^ ; 77 7 5 j 

2e/ + ec + 0(e/ - ec) 

where (p is the concentration of the crystal domains. This 
model also contains a single fitting parameter, cj). The 
mean field approximation is only accurate up to 0 ^ 0.2 
[17], but in the following we allow (j) to go all the way up 
to 1. 

Inserting the measured spectra of ec and e/ in Eq. 3 
does not produce a frequency shift. Consequently, this 
model cannot account for what we observe either. 

Instead, we propose to combine the two models such 
that a crystal layer is growing from the electrodes, while 
spherical crystallites are forming in the remaining liquid, 
see Fig. 4. This is modeled by combining Eq. (2) and 
Eq. (3) such that ei in Eq. (2) is given by the composite 
dielectric constant from Eq. (3). This model has two pa¬ 
rameters: the relative thickness of the crystal layer, ^siab, 
and the concentration of crystal spheres in the liquid, (j). 



EIG. 4: Cartoon of the combined heterogeneous and 
homogeneous crystallization. The growth is indicated 
by the arrows. 

Examples of fits to isothermal crystallization spectra 
for cell A and cell B are shown in Eig. 5(a) and (b), 
respectively. The fits are focused on the Debye process 
by only fitting to the points within the two dotted lines. 
Having two fitting parameters gives sufficient flexibility 
to account for both the decrease in relaxation strength 
as well as the change in peak frequency of the Debye 
peak. However, the combined model does not adequately 
explain the behaviour of the entire spectrum; it does not 
capture the behaviour of the alpha and beta relaxation 
processes during the crystallization or the broadening of 
the Debye process. Despite these limitations, the models 
ability to describe the behaviour of the Debye relaxation 
strength and peak position may be used to shed some 
light on the observed phenomena. 

The total crystallized volume fraction according to the 
model can be calculated as Xc = (1 — 0siab)0 + 0siab 
and the volume fraction taken up by the spheres alone 
as Xsphere = (1 “ ^siab)^- Figure 5(c) and (e) shows 
the volume fraction of the crystal layer and the spherical 
crystallites as well as the total crystal volume fraction. 
For both cells it seems that the onset of crystallization is 
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FIG. 5: The measurements shown in Fig. 2 analyzed 
using the Maxwell-Wagner (MW) approach described in 
section IIB. (a+b) show a selection of curves together 
with fits (magenta lines) for cell A (a) and cell B (b). 
The dashed vertical lines mark the frequency interval 
used for the fit. (c+e) show the volume fraction taken 
up by the spherical crystallites (circles) and the crystal 
slab (squares) as a function of time for cell A (c) and 
cell B (e). The solid line is the total crystal volume 
fraction, (d+f) Total crystal volume fraction from the 
MW fit (black line) as well as the normalized relaxation 
strength of the Debye process (in green) and alpha 
process (in red) as a function of time for cell A (d) and 
cell B (f). 


dominated by nucleation and growth of crystal spheres, 
and when a large fraction of the sample has crystallized 
(V r\j 80-90%), the slab growth takes over. However, the 
growth of crystal spheres starts earlier but proceeds at 
a slower rate in cell B compared to cell A. This differ¬ 
ence in crystallization behavior in the two cells that is 
reproduced at all temperatures (see Fig. 6). 

The total degree of crystallinity as calculated from the 
MW fits is plotted together with the normalized Debye 
and alpha relaxation strengths for cell A in Fig. 5 (d) 
and for cell B in Fig. 5(f). Again, the curves are clearly 
different for the two cells; in cell A Debye and alpha re¬ 
laxation strengths give almost identical curves that agree 
qualitatively with crystal fraction obtained in the MW 
fit, although the relaxation strength decreases faster than 
crystal fraction initially. In cell B, all three curves starts 


out in the same way but separates later in the process, 
where both Debye and alpha relaxation strengths overes¬ 
timates the degree of crystallinity. The MW crystallinity 
curve has a kink occurring approximately when there 
is a bump in the alpha relaxation strength, while this 
two-step behavior is not clearly seen in the Debye relax¬ 
ation strength. The behavior of crystallization process 
at 127.5 K and the MW analysis demonstrated in Fig. 5 
is general for all the studied temperatures as can be seen 
in Fig. 6. 

The proposed MW analysis qualitatively and quanti¬ 
tatively agrees with using the decrease of the alpha relax¬ 
ation strength rather than the Debye relaxation strength 
as a measure of the crystallinity of the sample. However, 
neither reflect the fact that sample does not crystallize 
fully. This is because both the Debye and the alpha pro¬ 
cess vanish during crystallization, while only the beta 
process survives. The analysis does not account for that, 
a fact that is already clear from the fits in Fig. 5(a+b). 


C. Avrami analysis 

Another independent - and more routinely used - way 
of evaluating the crystal growth is through the Johnson- 
Mehl-Avrami-Kolmogorov (JMAK) equation originally 
developed by Avrami [30-32]. In this analysis, the vol¬ 
ume fraction taken up by crystallites, Xc is expressed in 
terms of a growth rate constant k, an induction time to, 
and the Avrami parameter n as follows 

Xc{t) = 1 - exp [-k {t - to)"] . (4) 

The value of n depends on the crystal morphology and 
crystallization mechanism, but it is not straight forward 
to interpret the meaning of this parameter. Originally, 
it was a number between 1 and 4 such that n = d + 
r with d being the dimensionality of the growth and r 
being a number that represents the nucleation rate. For 
a constant nucleation rate r = 1 and r = 0 if nucleation 
stops when the crystallization starts. More recently n 
has been found to be a number between 1 and 7 [11, 15]. 

When dielectric spectroscopy is used to study crystal¬ 
lization the common practice is to assume that the alpha 
relaxation strength roughly corresponds to the degree of 
crystallinity, Xc, and then use that for the Avrami anal¬ 
ysis. Since the MW model also gives some support for 
the alpha rather than the Debye relaxation strength is 
expressing the degree of crystallinity in the sample, we 
will adopt this approach. 

One way to obtain an estimate of the parameter n 
is through the Avrami-Mehl-Johnson plot, which plots 
In [—ln(l — Xc(t))] versus ln(t). This procedure avoids 
fitting Eq. (4) to data and n is directly obtained as the 
slope of the curve. The Avrami-Mehl-Johnson plots for 
cell A and B is shown in Fig. 6. In both cases, we see a 
transition from a relatively high value n ~ 4 — 5 to a low 
value n ~ 1 — 2. This observations suggests a change from 
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FIG. 6: Results from MW and Avrami analyses for all temperatures, (a+b+d+e) Volume fraction of crystal spheres 
(a+d) and volume fraction of a slab (b+d) based on the proposed combination of Eqs. (2) and (3). The differences 
for the two measuring cells outlined in Fig. 5(c+e) are general; cell A (a+b) display a single peak in the “sphere 
growth” and the “slab growth” set in early, while in cell B (d+e) there is a clear double peak in the “sphere 
growth”, (c+f) show the Avrami-Mehl-Johnson plot which is a linearisation of the JMAK equation (Eq. (4) in the 
text). Eor both cell A (c) and B (f), the data clearly show a cross-over from one linear behavior with a steep slope 
4) at short times to another linear behavior with a smaller slope 1-2). But the transition from one to the 
other occurs at different crystallization degree for the two cells: for cell A the transition happens when the alpha 
relaxation strength has decayed more than 95% and for cell B when Ae(t)/Ae(0) ^40-50% (marked by the gray bars 

in the figure). 


higher dimensionality of growth to lower, which is con¬ 
sistent with the MW analysis suggesting a change from 
spherical to slab growth. Moreover, the data suggest that 
this transition happens earlier in the crystallization pro¬ 
cess for cell B compared to cell A, which could explain 
why the crystallization slows down and takes much longer 
in cell B. 


III. DISCUSSION 

Both the Debye process and the alpha process vanish 
during the crystallization, while the beta process sur¬ 
vives. Thus we confirm earlier findings that the crystal¬ 
lization process at temperatures near Tg stops before the 


sample is completely crystallized. Hedoux et al report 
signs of an aborted or frustrated crystallization process, 
signaled by a amorphous halo persisting in the x-ray spec¬ 
tra [21]. This slow and frustrated crystallization process 
has also been interpreted as a polyamorph transformation 
between two meta-stable liquid phases [24, 33-35]. Based 
on the dielectric spectra presented here it is perhaps dif¬ 
ficult to distinguish between the two scenarios, but the 
fact that the structural relaxation peaks disappear en¬ 
tirely combined with the emergence of Bragg peaks as 
documented in Ref. [21], point to a non-trivial crystal¬ 
lization process as the most obvious explanation for the 
observations. 

It is however interesting that the aborted crystalliza¬ 
tion is seen in the dielectric spectra as the survival of 






















the beta-process. If we envision the end product as a 
frustrated crystal, unable to tile space, then the liquid 
signal - in our case the beta relaxation - could originate 
from small pockets of liquids between crystal grains. This 
picture supports the idea of the beta relaxation being a 
local phenomenon, in favor of the “islands of mobility” 
suggested by Johari and co-workers [36, 37]. 

The mono-hydroxyl alcohols in general are interesting 
because of their anomalous (and usually intense) relax¬ 
ation process at frequencies lower than the structural al¬ 
pha relaxation - the so called Debye relaxation - which is 
believed to be due to supra molecular hydrogen bonded 
structures in the liquid [18, 38]. We observe that the De¬ 
bye process vanishes faster than the alpha during crys¬ 
tallization, and that the alpha intensity seems to give a 
better measure for the degree of crystallinity. Sanz et al 
[39] made similar observations for another monohydroxyl 
alcohol. They studied crystallization of isopropanol in 
real time by simultaneous dielectric spectroscopy and 
neutron diffraction measurements, and thus had a direct 
measure of the degree of crystallinity that could be cor¬ 
related with the relaxation strength of the Debye and 
alpha process. They observed that the Debye intensity 
dropped rapidly at the onset of crystallization, while the 
alpha intensity followed the crystallization. Their intu¬ 
itive and appealing interpretation was, that the breakage 
of the hydrogen-bonded network is a precursor of the 
crystallization, and that the molecules leaving the net¬ 
work did not immediately go into a crystalline structure. 
The MW polarization effects lends itself to a different 
- macroscopic - interpretation of the observations. Irre¬ 
spective of how the crystal growth is modeled in the MW 
framework, there cannot be proportionality between di¬ 
electric intensity and liquid fraction in the sample. The 
deviation from linearity depends on the specific model 
for the growth morphology and on the intensity of the 
process: the higher intensity, the stronger the deviation 
from linearity. Thus the MW analysis provides a sim¬ 
ple explanation for why the most intense process vanish 
before the less intense one. The MW analysis does at 
the same account for the observed frequency shifts of the 
relaxation processes. 

On the basis of MW fits, we suggested that the ob¬ 
served difference in crystallization behavior between the 
two cells could be rationalized by a transition from hav¬ 
ing primarily a growth of crystal spheres, a homogeneous 
nucleation and growth, to a growth of a crystal slab. 
This could be a slab growing from the electrodes, but 
could also be a certain point in the process where crys¬ 
tal grains percolate and effectively create a crystal layer 
in the liquid-crystal mixture. The difference between the 
cells would then be explained by a difference in the degree 
of crystallization when this transition takes place. The 
idea that the crystal growth changes from a higher dimen¬ 
sional growth to a low dimensional growth was supported 
by the Avrami-Mehl-Johnson analysis that also points 
to such a transition taking place at different crystalliza¬ 
tion degrees in the two sample cells. The overall validity 


of JMAK equation has been questioned, see e.g. Refs. 40 
and 41, and of course we need to be cautious when mak¬ 
ing conclusions, based on the MW analysis where the 
limits of applicability of the mean-field approximations 
was pushed. But since both types of analyses point with 
this picture, we believe that the proposed conception of a 
change in morphology of the crystal growth is consistent 
and sound. It remains to be shown how general this be¬ 
havior is. It would be interesting to apply this procedure 
to a simpler sample to study the influence of sample cell 
geometry on the crystallization process. 

Irrespective of the generality of the particular behavior 
found here, our study shows that one needs to be very 
cautious about making detailed microscopic interpreta¬ 
tions of the crystallization mechanisms based on dielec¬ 
tric spectroscopy alone, because MW polarization effects 
of the mixed phase requires knowledge about the crystal 
growth morphology. In addition, we have also shown that 
the crystallization is extremely sensitive to the specific 
sample environment. Thus, it would require extensive in¬ 
vestigations of different environments and perhaps even 
different probes to disentangle microscopic from macro¬ 
scopic effects. 


IV. CONCLUSION 

We have studied isothermal crystallization process in 
the deeply supercooled region of the mono-alcohol n- 
butanol in real time at 15 different temperatures using 
dielectric spectroscopy. Two different sample cells have 
been used to look for the effects of the sample environ¬ 
ment on the crystallization process. We found that the 
time evolution of the relaxation strengths differs for the 
two cells in a consistent and reproducible way for all tem¬ 
peratures. 

On the basis of the Maxwell-Wagner analysis, we sug¬ 
gest that the crystallization behavior can be explained by 
a transition from primarily growth of crystal spheres to 
growth of a crystal layer. The difference between the cells 
in this framework is the difference in when in the crys¬ 
tallization process this transition takes place. This pic¬ 
ture was supported by an Avrami-Mehl-Johnson analy¬ 
sis that also suggests a transition from higer dimensional 
growth to a lower one. 

The Maxwell-Wagner analysis can also account for the 
shift in peak frequency observed for the three processes 
during the course of crystallization, and thus a micro¬ 
scopic interpretation of the peak shift is not needed. 
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